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ABSTRACT 

At e + e~ and ep colliders, the large fluxes of almost on-shell photons accompany- 
ing the lepton beams lead to the photoproduction of jets. As the centre-of-mass 
energy is increased, regions of smaller x in the parton densities are explored and 
these are regions of high parton density. As a result, the probability for more 
than one hard partonic scattering occurring in a single 77 or jp collision can 
become significant. This effect has been simulated using an eikonal prescrip- 
tion combined with the HERWIG Monte Carlo program. The possible effects 
of multiple hard interactions on event shapes and jet cross sections have been 
studied in this framework at a range of energies relevant to HERA and LEPII. 
The results indicate that the effects could be significant. 
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1. Introduction 



For both protons and photons, QCD predicts a rapid increase in parton densities 
at low 1. In a naive treatment, this rise can lead to a corresponding (but ultimately 
unphysical) rise with increasing energy of perturbative QCD calculations of the jet 
contribution to the total cross section. However, the large number of small x partons 
contributing to jet production can mean that there is a significant probability for more 
than one hard scatter per jp or 77 interaction. The effects of multiple interactions 
can provide a mechanism for taming the rise in the QCD cross section^ in accord with 
unitarity. Clearly there are implications for the hadronic final state and in order to 
study these effects an eikonal model has been implemented within the hard process 
generation of HERWIG lift 

Multiple parton scattering, as shown in Figfj], is expected to affect jet rates. The 
average number of jets per event should be increased when partons from secondary 
hard scatters are of sufficiently high px to give jets in their own right. In addition, 
lower pt secondary scatters produce extra transverse energy in the event which con- 
tributes to the pedestal energy underneath other jets in the event. Thus multiple 
scattering can influence jet cross sections even when no parton from the secondary 
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Figure 1: An example of a multiple scattering in a 77 collision. 



scatters is itself of a high enough p T to produce an observable jet. By boosting the 
transverse energy of jets in this way, multiple scattering leads to an increase in jet 
cross sections for jets above a certain E^ 1 cut, even though the total cross section is 
reduced. 

In the model implemented in HERWIG, the number of hard scatters is given by 
the simple probability theory formula, 

^events = ^trials x event probability per trial. 
In the general case of a collision between particles (A and B) at an impact parameter 
b, the 'number of trials' is the product of the number densities of the partons in the 
overlapping regions of the parent particles, i.e. 

#trials = J d 2 b' n i/A (b', x A )n j/B (b' - b, x B ) = A(b) n i/A (x A )n j/B (x B ), (1) 

where ni/ A (x A ) is the number density of partons of type i within particle A and 
carrying a light cone momentum fraction, x A . We have assumed that the x and b 
dependences of the parton densities factorize (we call A{b) the area overlap function). 
The '^events' is the multiplicity of final state jet pairs produced per interaction. The 
'event probability per trial' is daij/(PX 1 P B 1 crH(s)) where a B (s) is the cross section 
for the production of one or more jet pairs and d&ij is the parton-parton subprocess 
cross section. P A is the probability that the particle A actually interacts as a hadron. 
For hadrons, the factor is unity, whilst for photons it is proportional to a cm . In fact 
we take the p-dominance form P 1 = 4na eTa / p p . 

The b dependence of the parton density is assumed to be given by the electro- 
magnetic form factor of the parent particle. By normalizing the integral of the area 
overlap function to unity we can identify 

A (b)d%^fP- = X (b, s)d% (2) 
r A f B 



with the mean number of hard scatters occuring in the element d 2 b of impact pa- 
rameter space (per hadronic interaction). This defines the 'eikonal' function, x{b, s), 
in terms of <J^ C , the total inclusive cross section for the production of jets (with 
Pt > PTmin) which is given by, 



ps/i pi pi 

a%°(s) = dp 2 T dx A / dx B J2fi/A(x A )fj/B(x B 



(3) 

where the parton distribution functions, fi/A( x A), are related to the number densities, 
ni/A(xA), by fi/A(xA) — PA n i/A( x A)- This inclusive cross section includes the mean 
multiplicity of hard scatters per event, i.e. tf"# c (s) = (n H )aH{s). This ensures that 
<7ff(s) cannot be larger than the total cross section, whereas, in principle, a l £ c (s) can 
be. Assuming the successive scatters to be uncorrelated, Poisson statistics then give 
the probability for m (and only m) jet pairs to be produced: 

Pm= ix(b,s)r (4) 

ml 

The HERWIG Monte Carlo generates the required number of hard scatters and 
the associated initial and final state parton showering. The outgoing partons and 
remnant jets are then fragmented to the hadronic final state. Two modifications are 
made to the simple eikonal model in the implementation. Firstly, energy conservation 
is imposed, i.e. after the backward evolution of all the hard scatters in an event, the 
energy remaining in the hadronic remnants must be greater than zero. Secondly, if 
during the backward evolution of the first scatter, the splitting qq <— 7 is arrived at 
before the evolution cut off scale, the event is classified as an 'anomalous' event and 
no multiple interactions are allowed (this is a conservative approach). 

2. Multiple Interactions at LEPII 

The CM energy at LEPII is sufficiently high that the model predicts a significant 
proportion of multiple scatter events. This can be seen in table [1], where we show the 
expected mean number of hard scatters per event at a variety of colliders. In each 
case, we defined a hard scatter to have px > PTmin and used the MRS D_ i proton and 
GRvi photon parton densities. Note that we have generated only 'twice resolved' 77 
events. For these events, since the cross section falls rapidly with increasing pj>, we 
expect that most of the multiple interactions will have px ~ PTmin and so their main 
effect will be to increase the mean number of observed jets by boosting the underlying 
Et in the event. Events which contain multiple interactions with high enough px to 
be observed as jets in their own right are relatively rare, but their observation (they 
appear as pairs of back-to-back jets) would provide striking evidence in support of 
the existence of multiple interactions. 
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p 

1 Train 




Mean no. scatters 


LEPII 


2.0 


180 


1.123 


LEP 


2.0 


90 


1.084 


TRISTAN 


2.0 


60 


1.080 


HERA 


3.0 


296 


1.04 



Table 1: The mean number of hard interactions per event. 



Jet finding was performed on the hadronic final state using a cone algorithm! with 
cone radius R = 1. Jets have Ex > 3 GeV and pseudorapidity —3 < "\f et < 3. The 
inclusive jet cross section, as shown in Fig.^a, is only affected by the jet boosting, i.e. 
for a PTmin = 2 GeV, scatters from 2-3 GeV are boosted over the E T e cut. The effect 
of multiple interactions is most significant in higher jet multiplicity rates, as can be 
seen in Fig.|]b and c, where the dijet cross section is shown. 

A variable found to be useful at HERA i for distinguishing between direct and 
resolved photon events is x° bs . It is defined by 

and the sum is over the two highest Ex jets in the event. In leading order, this is 
the fraction of the photon energy which enters the hard scatter, i.e. direct events 
are peaked at x° bs = 1 whilst resolved events have x° bs < 1. The same variable can 
be used in 77 collisions where twice resolved events populate the lower x° bs region. 
In Fig.^Jd, the cross section is plotted as a function of x° bs . It can be seen that the 
effect of multiple scattering is greater in the low x° 6s region. Multiple scatterings are 
more likely to occur here as it is in this region that the higher parton densities occur; 
also the energy conservation constraint is less restrictive. We have not included once 
resolved and direct 77 events which will populate the high x° bs region and should be 
relatively unaffected by multiple interaction. 



3. Multiple Interactions at HERA 

For HERA, the minimum transverse momentum of a hard scatter is set to prmin = 
3 GeV. A cut on the 7p CM energy was made, i.e. 114 GeV < yfs < 265 GeV, similar 
to those usually made by the experiments. For these choices the mean number of hard 
scatters per event was found to be 1.04. 

Jet finding was again performed with cone radius R = 1. Jets have Ex > 6 GeV 
and pseudorapidity —2 < rf et < 2. The additional Et can be seen directly in the jet 
profile, Fig.[|a, where the E T in the jets is plotted against rj relative to the jet axis. The 
pedestal energy is increased with the inclusion of multiple interactions. The x° 6s cut 
isolates the resolved interactions, and as expected multiple interactions have no effect 




Figure 2: a) E^ 1 distribution for inclusive jets b) E^ 1 distribution for dijets c) rf et 
distribution for dijets d) ^^distribution for dijets, at LEPII energies, with multiple 
scattering (solid line) and with multiple scattering turned off (broken line). 

in the direct case, Fig.^b. The extra E? enhances the inclusive jet rate around 
Yjjet _ -j^ an( ^ an i ncreasec l sensitivity to multiple scattering can be seen in the rf et and 
Ej? distributions in dijet events, which are shown in ref.4. An enhancement at high 
rf et and low E^ is seen. The x° bs distribution, Fig.|3]c, has the direct contribution 
generated using HERWIG included, hence the rise at x°^ s ~ 1. The inclusion of 
multiple scattering has a significant effect on the lower a;° bs region, as in the 77 case. 

4. Conclusions 

The effect on the hadronic final state of multiple parton scattering in 77 and 7p 
interactions has been simulated by interfacing an eikonal model of multiple parton 
interactions with HERWIG. Multiple scattering must occur at some CM energy, in 




Figure 3: The Et jet profile in 77, for a) x° bs < 0.75 and b) x° bs > 0.75. c) 
^^distribution for dijets, at HERA energies, with direct contribution. Including 
multiple scattering (solid line) and with multiple scattering turned off (broken line). 

order that unitarity is not violated. Our model indicates that the effect of multiple 
scattering is significant at both HERA and LEPII energies. For reasonable experimen- 
tal cuts, the inclusion of multiple scattering leads to significant changes in inclusive 
and dijet cross sections which should be understood before attempting to unfold to 
parton distribution functions. 
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